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ABSTRACT 

This paper is concerned with the creation and subsequent motion of singularities of 
solution to classical Rayleigh-Taylor flow (two dimensional inviscid, incompressible fluid over 
a vacuum). For a specific set of initial conditions, we give analytical evidence to suggest the 
instantaneous formation of one or more singularity(ies) at specific point(s) in the unphysical 
plane, whose locations depend sensitively to small changes in initial conditions in the physical 
domain. One-half power singularities are created in accordance with an earlier conjecture; 
however, depending on initial conditions, other forms of singularities are also possible. 

For a specific initial condition, we follow a numerical procedure in the unphysical plane 
to compute the motion of a one-half singularity. This computation confirms our previous 
conjecture that the approach of a one-half singularity towards the physical domain corre- 
sponds to the development of a spike at the physical interface. Under some assumptions 
that appear to be consistent with numerical calculations, we present analytical evidence to 
suggest that a singularity of the one-half type cannot impinge the physical domain in finite 
time. 


‘Research was supported by the National Aeronautics and Space Administration under NASA Contract 
No. NAS1-18605 while the author was in residence at the Institute for Computer Applications in Science 
and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. INTRODUCTION 

The motion of the interface of a heavy fluid initially resting on top of a lighter fluid 
(Rayleigh-Taylor flow) is a very basic but important problem in fluid dynamics and has 
been the subject of intensive research over a long period of time. Recent interest in 
the Rayleigh-Taylor instability stems from its disruptive presence in inertial confinement 
devices (See Verdon et al (1982) for instance ). Emmons, Chang k Watson (1959) studied 
the interfacial features experimentally with initially sinusoidal disturbances and found 
that for large times, a pattern containing downward spikes and upward moving bubbles 
forms. In other experiments (Read, 1984) with Atwood ratio close to one (i.e. density ratio 
between lighter and heavier fluid close to zero), a variety of bubbles and spikes is formed for 
random initial condition. It is clear that there is. a significant interaction between bubbles 
(or spikes) so that the congregate motion is rather different from a regular pattern. 

In the idealized limit of two-dimensional in vis cid, incompressible fluid over a vacuum, 
direct numerical calculations by Baker, Meiron k Orszag (1982) have shown that an initial 
sinusoidal perturbation of the interface leads to an upward moving bubble and a downward 
moving spike in each period of the disturbance. The shape of the upward moving bubble 
agrees with the steady bubble solutions of Davies k Taylor (1950), while a downward 
moving spike accelerates with free fall. For initial conditions containing more than one non- 
identical undulation per period, the demand for appropriate resolution makes it difficult 
to continue calculations (along the lines of Baker, Meiron k Orszag (1982)) for sufficiently 
long-times to identify the effective acceleration of the bubble-tip-envelope observed in the 
Read experiment (1984). 

Given the physical importance of its dynamics, the Rayleigh-Taylor problem has also 
been studied from a more practical perspective since direct numerical calculations based 
on the fluid-dynamical equations appear to be impractical even for the simplest of the 
Rayleigh-Taylor flows. Model equations have been developed (see Gardner et al (1988) 
and Sharp (1984) and references there in) to study the interaction of multiple bubbles and 
spikes. Typically these include parameters that in some cases can be computed rationally 
by appeal to physical dynamics. Some of the more recent models (Gardner et al (1988)) 
have been developed in great generality without any of the restrictive assumptions of an in- 
viscid incompressible two dimensional flow. While these studies have been quite important 
in furthering physical understanding of bubble competition and merger processes, we are 
unaware of any direct derivation of model equations from the fundamental fluid equations. 
While this paper does not address this problem either, we hope the approach in this paper 
will eventually bridge the gap between the direct numerical simulation of fluid equations 
and model studies, at least in simple cases. 
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Here, we explore the dynamics of singularities in the classical Rayleigh- Taylor problem 
without resort to any localized approximation (Siegel (1989), Baker, Caflisch &; Siegel 
(1992)). At time t, consider the conformal map z(£,t) that maps the interior of a cut unit 
circle in the ( plane (Fig. 1) into a periodic strip in the physical domain (Fig. 2) such 
that the origin coincides with z = —i oo. The unit circular boundary then corresponds 
to the free boundary. We will be concerned with the formation and subsequent motion of 
singularities of z(£, t) and the complex velocity potential IF(£,t) in the unphysical domain 
|C | > 1. We have several long range goals in furthering such an understanding. 

First is the possibility that singularities can be analytically subtracted out in a basis 
representation of / and W making them amenable to direct numerical calculations for a 
long time. Second is that the bubble and spike interaction can be understood through 
the interaction of multiple singularities in the unphysical complex plane. As shown in this 
paper, a one-half singularity in |£| > 1 approaching |<C| — 1 corresponds to a continually 
developing spike at the physical interface. The portion of the unit circle |£| = 1 between 
any two approaching singularities contains the image of the bubble boundary in the ( 
plane. Thus, pairwise singularity merging corresponds to a bubble getting smaller, while 
its neighbor becomes larger a well known process in the Rayleigh- Taylor problem. Third, 
it may be possible to reach general qualitative and quantitative conclusions about the 
relation of long time bubble dynamics to the specifics of initial conditions in the complex 
unphysical plane, which is related to the physical initial condition in an ill-posed way. This 
may allow one to construct an appropriate statistical model of bubble interaction in terms 
of the statistical features of the initial conditions in the unphysical domain. 

In a previous paper (Tanveer, 1991a), the analytically continued equations for the 
two-dimensional Rayleigh- Taylor and water wave problems were derived in the unphysical 
domain |£| > 1. For steady water waves, analytical and numerical calculations were 

carried out to establish the relation of water wave crests to one-half singularities of z(£, t). 
However, for the Rayleigh- Taylor or the unsteady water wave problems, no concrete results 
were obtained except to note that certain one-half singularities of z(£, t) and IF(C> t) were 
consistent with these equations. It was also noted that in the limit of a one-half singularity 
approaching the physical domain, the analytically continued acceleration at a one-half 
singularity is the free-fall under gravity, similar to that which is observed for a spike. 
Based on this, it was conjectured that a one-half singularity approaching the physical 
domain corresponds to a spike developing at the physical interface. The work presented 
here is is a natural continuation of our previous (Tanveer, 1991a) work. This paper is 
organized as follows. 

In Section 2, we present the analytically continued equations in the unphysical plane 
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|C| > 1 that has been derived previously (Tanveer, 1991a). The equations are presented 
in several alternate forms, some more convenient for asymptotic analysis, while others for 
numerical computation. 

In Section 3, we show that under some assumptions on the single-valuedness of ( as 
a function of a characteristic variable in some region of the characteristic space, the only 
possible singularity ( s (t) of z(( y t) and in ( is of a “fold” type where each of dL and 

7^- are analytic functions of the variable (( - C*(0) 1/2 or (C - C*(0) 1/3 or (C “ C*(*)) 1/4 > 
etc. In the case when the fold singularity is of the one-half type, we relate the coefficients 
in an expansion in [(“ — Cs(f)] 1 ^ 2 to the solution in the characteristic plane. 

In Section 4, we address the question of singularity formation-how does a singularity of 
z (Ci 0 an d W (C, t ) form in the unphysica.1 domain when there is none initially? We consider 
several classes of initial conditions for which z(£,0) and W(£,0) are analytic everywhere 
in the finite £ plane outside the unit circle. We give analytical evidence to suggest that 
singularities can form instantaneously at a point in the £ plane where Z£((\0) = 0. This 
result is very similar to results obtained in similar situations for other fluid flows such as the 
Hele-Shaw flow with surface tension (Tanveer (1991b)) or the Kelvin-Helmholtz problem. 
In the latter case, recent work of Cowley et al (1992) has shown that the Moore singularity 
(Moore (1979,1985), supported by numerical computations of Krasny (1986), Shelly (1992) 
and rigorously analyzed by Caflisch & Orlenna (1988), actually forms instantaneously at 
some point in the complex circulation variable. We also find that in our problem, for 
some initial conditions, a singularity moves in instantaneously from infinity to the finite ( 
plane in the sense that for any fixed t > 0, the singularity is at a finite ( point; yet as 
^ > 0"*") this location recedes to infinity. Our calculations suggests that for certain set 

of initial conditions, only one-half singularities can be created; however, there exists other 
initial conditions for which singularities of a more complicated form involving logarithm 
can occur. In this case, the assumption on single valuedness of £(£,<) that is assumed in 
the analysis of Section 3 is violated. 

In Section 5, we employ a numerical procedure to track the motion of a one-half 
singularity that is created at the initial instant of time. We compute not only the location, 
but also a few coefficients of the one-half power expansion. So far, numerical computation 
has been performed for a very special initial condition. Nonetheless, the result confirms our 
previous conjecture that the approach of a one-half singularity corresponds to a continually 
developing spike at the physical interface. 

In Section 6, we address the question if an approaching one-half singularity of z(£,<) 
and W((,t) of the type computed in Section 5 can actually impinge the physical domain 
boundary |C| = 1 in finite time. With certain assumptions that appear to be consistent 
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with the numerical calculations in Section 5, our analytical evidence suggests that an 
isolated one-half singularity cannot impinge the real domain in finite time. However, this 
leaves open the possibility of different kinds of singularity or multiple one-half singularities 
coalescing at |(| = 1. This result has a bearing on the work of other researchers. Siegel 
(1989) and Baker, Caflisc.h & Siegel (1992) have studied exact complex travelling wave 
solutions to a localized simplification of the Rayleigh-Taylor equations for arbitrary Atwood 
ratio. For unit Atwood ratio (the case studied here), Baker, Caflisch & Orlenna (1992) 
found a class of travelling wave solutions with one-half singularities, each of which moves 
at a constant speed. However, based on a spectrum fit of the numerically computed results 
for the full Rayleigh-Taylor problem in the physical domain, they detect a definite slow 
down of such singularities at unit Atwood ratio which is at variance with the solution to 
the localized approximation. However, it remains unclear from their work if the slowdown 
was sufficiently significant to avoid a finite time singularity in the real domain. 

Since the analytic continuation of z((, 0) from |(| < 1 to the unphysical plane |<| > 1 
is an ill-posed procedure, i.e. arbitrary small deviations of z(C 5 0) in |C| < 1 can affect the 
location of its zeros of 2^(C, 0) in |C| > 1 5 it follows that the precise location and number 
of such singularities created at initial time will be highly sensitive to initial conditions in 
the physical domain. Our results on the correspondence of singularities with spikes at later 
times can explain the observed random nature of bubble-spike interaction in the long-term 
behavior of the physical interface. In our discussion in Section 7, we make a plausibility 
argument on how singularity interactions can explain bubble competition. 

2. MATHEMATICAL EQUATIONS 

The conformal map from the cut unit £ circle (Fig. 1) into a periodic, strip in the 
physical domain as shown in Fig. 2 ( z = x + iy ) can be decomposed into 

z((,t) = 2tt + ilnC + */(C»0 ( 21 ) 

where /(£, t) is oblivious to the branch cut and therefore possesses a convergent power 
series representation for |£| < 1 

OC 

/« ,i) = E«.(')c < 2 - 2 ) 

n=0 

Here we have assumed, without any loss of generality, that the period in the z plane is 27r 
and the acceleration due to gravity is unity and is directed upwards (along the positive y 
axis). For analytic shape, the convergence of (2) occurs up to |Cl = 1 ■ Similarly, there 
exists a power series representation for the complex velocity potential 

OO 

W((,t) = MOC" • ( 2 - 3 ) 

n-0 


4 


We will assume that the initial conditions are symmetric, so that a n and b n are initially 
real. From the equations, it is clear that these symmetries are preserved for all later 
times. This assumption is only made for simplicity and generalizations for nonsymmetric 
disturbances are possible. This means that on the real ( -axis in the interval (-1,1), 

Irrt. / = 0 (2-4) 

holds for / and the complex velocity potential W satisfies 

ImW = 0. (2.5) 


The kinematic condition on the free boundary can be expressed as 

^ In p(x,y,t) = 0 (2.6) 

on p(x,y } t) = 1 , where £ = p e lu , with v real. In this representation, In p , v and 
t can be thought of as three dependent variables depending on x , y and t . Switching 
the role of dependent and independent variables, the kinematic condition implies that 


Re [ C W c - C* z l *t] = 0 


(2.7) 


where the symbol * here and in what follows stands for complex conjugation. Plugging in 
the representation for z from (2.1) on |£| = 1 , we find that (2.7) is equivalent to 


Re [ 


ft 

1 + < k 


Re CW f 

|i + C/cl 2 


( 2 . 8 ) 


on ( = e tv for v in the interval [0, 2 7r] . The analytic continuation of this for |£| > 1 
(see Tanveer, 1991a for details) is 


u e w< + i tv c (i/c,q _ 

i + C /c (1 + C /<) (1 + ^/c(i/C>0) 

where I 2 can be written as either of the following two expressions 


MC ,0 - 


1 + 


(2.9) 


( 2 . 10 ) 


W,i) 


1 

4 wi 


/ 

•'1C' 1=1 


dC 

c 


rc + c'l C'WdC'J) + jr WcCjr.Q 

.C'-cJ [1 + C + ? 


( 2 . 11 ) 
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The Bernoulli’s condition on the free surface for this time dependent problem can be 
written as 

(2 ' 12) 

on the unit circle ( = e* " . The analytic continuation of this outside the unit circle (see 
Tanveer, 1991a for details) is 


W t - 


C w< ft 


f + 


W c (l/C,f) 


i + ( /< J (i + C /<) (i + ^ + /<( i/CiO) 

where / j can be written in either of the following two representations 

WWM) w<(;M 


(2.13) 


(2.14) 


_L 4 


+ 
1 

47TZ / 

/•/ 

IC' =1 

Lc-cJ 


(1 +C /<«',<)] (1 + f7 /<(£,<)] 

Equations (2.7) and (2.13) can be written in a more convenient form by defining 


(2.15) 


y i = 


V2 = 


i + C/c 

l 

TTOc' 


In that case, we get (Tanveer (1991a)) 

y lt — (i? 3 + i?2j/i)7/i c = Cyiy2-Ri< - (i + C-fi<)y2 + i 


(2.16) 

(2.17) 

(2.18) 


V2 t - (#3 + R2y\)y2 i = -R2V2V\ i + ^-y2-^-yl -C^yf- C(-4)<yiy2- #3<y2 H-^yf 

(2.19) 

(2.20) 


where 


Un^t) 1 

1 + ? <’ 

= 1 + Twpj = c !/2< ^’ i) 


ft = </ 2 

R = 1* + . < , „ = h + ».(^,«). 


1 + < /<(<,<) 


( 2 . 21 ) 

(2.22) 

(2.23) 
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Notice that 7i«,t) and /.(c. /) given by (2.15) and (2.10) can also be written as 



i 

+ 

y r> 

k' =i s 

Lc-cJ 


h(C ,t) = 


2ni 


/ f 

J IC'l=i 


C + C' 

Lc'-cJ 


|yi (C'.Of 


Re [y*i(C,t)y2(C,t)} 


(2.24) 

(2.25) 


since from symmetry properties (2.4) and (2.5) and the relations (2.16) and (2.17), y*(£' ,t.) - 
yi(l/£',i) and y£(£',*) = y 2 ( l/£',<) on IC'I = 1- By using (2.18) and (2.19), equation (2.19) 
can be replaced by a relatively more compact equation 


1 


Xu 2. 

t 


(i?3 + R 2 y \ ) 


Cy2 J 


—Ric . 


(2.26) 


By introducing an appropriate characteristic variable £ such that £ = C(£>0 and 
£(£,0) = £ and defining 


yi(£,<) = yi(C(£,0,<) 


hits) 


yXXXXhXl 

C<(£,0 


(2.27) 

(2.28) 


one finds that equations (2.18) and (2.19) is equivalent to the following set of equations 
for £(£,<), yi(£,0 and y 2 (£,0 (Tanveer (1991a)) 


0 = — i? 3 - i? 2 yi (2-29) 

yi t = (.inh Q R\ ( — (i +C^i c )y2 £$ + 1 (2.30) 

y 2 t = -^-£2 - -^r £2 C* - C -&ic C* £2 + ^3 ( £$ 2/2 + -jriiih- (2.31) 

Alternatively, from (2.29)-(2.31), we can derive 

-4-1 = -JMC(£,*),0C* (2-32) 

LCy2 J ( 

'§]=-[! + C *<(«{, 0.01 f + jj-. (2.33) 


Equations (2.29), (2.32) and (2.33) will form the basis of the numerical calculation de- 
scribed in Section 4. 
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3. PROPERTIES OF A CLASS OF SOLUTIONS 


Consider initial conditions for which each of ^((,0) and TF^(CjO) are analytic every- 
where in \(\ > 1 except possibly at oo. Also, assume that there is some open region 7Z 
in the £ plane in |£| > 1 so that the image of R under £(£,£) U P to certain time T is 

contained in \(\ > 1 with no point mapping to ( = oo. Further, we require that 
is single valued in 7Z. These set of assumptions will he referred to later as Assumption A. 
Our results in the next section suggests that Assumption A can only be valid for arbitrary 
K in |£| > 1 for some class of initial conditions. 

Consider an arbitrary closed contour C within 7Z where the Assumption A above is 
valid up to some time T > 0. We now derive some analyticity properties of the solution 
in this region up to time T. 

It follows from (2.32) that 


d 

dt 



= 0. 


Integration of (2.33) implies that 


d 

dt 



dt 


Vi 

C V2 



(3.1) 


(3.2) 


Initially, £(£,0) = t an d J/i an d $2 coincide with j/i and 3 / 2 • Thus since z^(^,0) and 
W^(4, 0) are each analytic (Assumption A), from (2.16) and (2.17), it follows that each of 
^ ^ and £<* is an analytic function of 4 for |£| > 1 (except possibly at 00 ). Thus, 

C (4 o)k (4 0) an< ^ C(/o)| 2(4 0) are anal y tic functions of £ for |£| > 1 (except possibly at 00). 
Thus, from Cauchy’s theorem 



1 

C(£,0) fc(£, o) 



yi(eo) 
aeo) w(e,o) 


From (3.1) and (3.2), it follows that 



1 


0 

0 . 

0 


(3.3) 

(3.4) 


(3.5) 



MU) 

tit,*) Mtu 


0 . 


(3.6) 


Then from Moerara’s theorem, ancl cl/bldh) ana ly tic in t in K. Then, 

yi(t,t) and y 2 can only have pole singularities of the same order. However, from (2.20), a 
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pole in yj is compatible only with a logarithmic or worse singularity of £(£>*)» which 
therefore violates Assumption A. Thus, under Assumption A, we find that each of y x and 
y 2 will be an analytic function of £ in the region R. 

However, despite the analyticity of y\ and y 2 as a function of £ under Assumption A, 
each of y, and y 2 (and therefore / and W) can have singularities in the image of R in 
the £ plane, as we shall now see. A singularity appears whenever there is failure of local 
inversion of the relation £ = <(£,*) into £ = £(£,<)■ Such a singularity will be referred to 


as a fold singularity as will occur a point £ 0 (<) where 

<{(?«(<),<) = °- < 3J) 

For a fold of the simplest kind, 

^ 0. (3-8) 

It is clear that if we define ( s (t) = £(&>(*)> near £ = 

c = ao + 5C«(eo(<),oie-MO] 2 + - 

C* = C«(£o(OiOl£ - &>(*)] + (3.io) 

y a (£,t) = ii(t) + A 2 (0(£ - 6(0) + ^3_(0(f - Co(0) 2 + O (£ _ 6>(0) 2 t 3 - 11 ) 

y 2 (£,t) = B 2 (t) + 5 3 (£-£o(<)) + 0(£-£o(0) 2 - ( 3 - 12 ) 

Then, it is clear from (2.27), (2.28), (3.9)-(3.12) that 

yi (C,0 = A,(0 + a 2 (0(C-C.(0) 1/2 + ^(C-C.) + °(C - (sf /2 ~ (3.13) 

and 

y 2 (C,t) = B 2 (t) (C-C.(*)) 1/2 + 5 3 (<XC-Cs) + o((C-C) 3/2 )- (3- 14 ) 

where 

A,(t) = ^i(0 ( 3 - 15) 

Mt) = \] MUtM Mt) (3 ' 16) 

B 2 (t) = yfe ■ ( 3 ‘ 17 ) 


Similar expressions can be found for the other coefficients. Going through the inversion 
process carefully, it is clear that the analyticity of y\ and y 2 at £ = £o(t) implies that each 
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of yi(C,0 and y 2 ((,t) is analytic in [( - C*(f)] 1/2 at ( = ( s (t). From the definition of y, 
and y 2 in (2.16) and (2.17) and that of f in (2.1), it follows that under Assumption A and 
condition (3.8), each of and 7F(£,t) will be analytic in [£ — Ca(t)] 1 ^ 2 at ( = (s(t)- 

The speed of such a singularity («, can be found by noticing that due to (3.7), 
dt C(£o(0> 0 Ct(£o(0> s ' nce C<(eo(i). 0 = 0 for a fold singularity. Therefore from (2.29), 

Cs = - R 3 ((s(t),t ) - R 2 (( a (t),t) y x (( 3 (t),t). (3.18) 

Alternatively, on substituting (3.13) and (3.14) directly into (2.18) and (2.19), we find 
that the coefficients of the most singular terms as ( — > ( s (t) is given by 

4 d > 

A 2 j~(s + A 2 R 3o + R 2g A 2 A\ = 0 (3.19) 

and 

d 

B *J t C. + (^3 0 + R^)B 2 = 0 (3.20) 

where subscript 0 refers to the evaluation of those quantities at ( = ( 3 (t) . Note that 
each of equations (3.19) and (3.20) are consistent with (3.18). Equating progressively less 
singular terms in ( - ( s (t) obtained by substituting (3.13) and (3.14) into (2.18), we get a 
set of relations, the first two of which is 


= 1 + \ A l ^(C-CO.O (3.21) 

^2 = 2 [-^3 T R 2( ;(£ s (t), t) A\ + t)Az\A 2 

+ R 2 (Ut),t)A 2 A 3 + C s AiB 2 R H (Cs(t),t) - (l + C./i<(C,(0»*))*2. (3.22) 

Similarly on substituting (3.13) and (3.14) into (2.26) and equating different powers of 

(C Cs) * we find that the leading order equation is just (3.20). At the next order, we 
find 


d 

1 

1 

dt 

[c,b 2 j 

2 


CsB 2 


(JMC.(0,0 + ^2c(C.(0.0^i + R 2 (Cs(t)A) A 3 ) 


R2((s{t ), t)A 2 


B, 


BUs 


(3.23) 


It is clear from this expansion that once A \ , A 2 , ( s (t ) B 2 (t ) and 53 (f) are known as 
a function of time, one can calculate the coefficients of the higher order terms (say A 3 
and B :i ) in terms of the already known lower order terms, provided global terms such as 
f^3(Cs(t),f) and 1 1 £(£*(<),<) are known. This is an important observation as it gives us, at 
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least in principle, a method of calculating as many coefficients in the one-half expansion 
in (3.13) and (3.14) as we want, provided the first few terms are known. In Section 5, we 
numerically calculate the first few coefficients and also show how the global terms can be 
calculated. 

Using a similar procedure as above, it is clear that if £o(t) is a double zero of Q 
and Assumption A holds, the expansion for and j/ 2 ( C ? ^ ) will contain powers of 

(( — ( 3 (f)) 1 / 3 . Generally for a zero of Q of n-th degree, one can expect a series in 
[£ — C«(t)] 1/,< ' n+1 ^ f° r each of yi and x) 2 - 


4. SINGULARITY FORMATION IN THE RAYLEIGH-TAYLOR PROB- 
LEM 

In this section, for certain initially analytic 2 (C) 0) (hence f((, 0)) and W(£,0), we give 
analytical evidence to suggest the instantaneous formation of one or more singularities in 
f(Cit) and TF(£,t) in the unphysical domain |£| > 1 at certain points that depend on 

properties of z(£, 0 ) and W(C,0) in |C| > 1 - 

First, note that, for 0 < t << 1, we can try a regular perturbation expansion in t. 
It is clear from (2.18) and (2.26) that 


yi «, t) = yi (C, 0 ) + 1 [{ J? 3 ( c, o) + fl 2 «, o)y, (C, o)} y u (C, 0) + Cyi (C, o)y 2 (C, o)i?i (C, 0 ) + 1 


1 


Cyz Cy2(C,o) 


t 


-[i+ai t (c,o)]ift(c.o)i] +o(( 2 ) 

2^3{C» 0) + H2(C,0)yi(C,0) 1 n /> 




+ 0{t 2 ). 


(4.1) 

(4.2) 


Clearly, since each of R\ (£, t ), I? 2 (C) 0> -^(C? 0 and (C> 0 involve y\ and j /2 in the physical 

domain £ < 1 (where each of them can easily be seen to have a regular perturbation 

series in i), one will have a regular perturbation expansion in powers of t, of these global 
quantities as well. On substituting into (2.18) and (2.26), the coefficients of all powers of 
t can be determined, at least in principle. 

This regular perturbation series for t << 1 becomes disordered when the later terms 
of the perturbation expansion is more singular than the previous terms, which can occur 
at some finite (0 or as ( — > 00 . Here, we will only restrict to breakdown in the power 
series due to one or more of the following conditions: 

(a) = 0 at C = Co with z cc (Co, 0 ) / 0, W<(Co, 0) 7 t 0 

(b) = 0 at ( = (0 with 2 <^((o, 0 ) ^ 0 , FF(C, 0 ) = 0 for all (. 

(c) As £ — » 00 , z^((,0) ~ & 2 C m_1 f° r positive m and W^(^,0) ~ b\ £( n + m_1 ) for n > 1 . 
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4a. Birth of a pair of one-half singularities 

Here we assume that condition (a) holds at some point Co- In this case, from the 
definition of / in (2.1) and y\ and y 2 in (2.16) and (2.17), it follows that the asymptotics 
°f yi (C> 0) and 2/2 (C? 0) as ( — > Co is given by 


yi(C,o) ~ 


1/2 (C,o) ~ 


f»2«(Co,0)(C - Co) 

(4.3) 

w'cCCo.o) 

(4.4) 

Co 2 c<(Co. 0)(c — Co) 


It is clear that the regular perturbation expansion in (4.1) and (4.2) get disordered as 
C * Co f°r any t > 0. Indeed the coefficient of t in (4.1) and (4.2) become the same 
order as the leading order terms when ( — Co = 0(t x / 2 ). This suggests the choice of 


_ C 1 — C/Co) 
1 1 1/2 


(4.5) 


as the inner variable. In the limit t -> 0+ with r/ = 0(1), one then finds from (2.18) 
and (2.26) that the solution is given in the similarity form 


yi = a -1 *“ 1/2 rift) 


(4.6) 


where 


V2=t- 1 / 2 Yifo) 

_ -^2 (Co ) Q) 

Co 


with a assumed nonzero. With 0(t 1 / 2 ) error, each of Y\ and Y 2 satisfy 


(4.7) 

(4.8) 


y 1 + v y;- 2 y 1 y; = 0 (4.9) 

y 2 + r ? y 2 ' + 2I2 y; - 2 r, k; = 0 . (4. 10) 

The asymptotic matching condition that (4.6) and (4.7) match to the initial condition that 
behave like (4.3) and (4.4) near the singular point implies that we must require that for 
large \i]\ 

, x ab 

Y\(j]) ~ — (4.11) 

Y 2 (r}) ~ C - (4.12) 
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where 


and 


b= - 


W(« 0 ,0) 


Co/«(Co,0) + / c (Co,0) 

1 


c = 


Co 2 /«(Co,0) + Co/c(C«,0)' 

Equation (4.9) and (4.10) with conditions (4.11) and (4.12) can be solved exactly 


(4.13) 

(4.14) 


r,w = l-J£-*a 


n(i) 


4b 2 , i 2 


\/,, 2 — 4 ba[r) — \J r/ 1 - 4 ba] 


(4.15) 

(4.16) 


The above solutions (4.15), (4.16) show that two square root singularities at 7 / — 
±\/4 6a, i.e. at ( = Co T Co \/4 bat for 0 < t << 1 for each of and y 2 eac h °f 

which is consistent with the expansion (3.13) and (3.14) for sufficiently small |Co TCo V^bat\. 
Going to the definition (4.6) and (4.7), it is clear that for t << 1, each of A 2 (t), 

i?2(0 are directly determinable from (4.15) and (4.16) at each singularity. 

4b. Singularities involving logarithms 

Here, we assume condition (b) at some point Co- With this initial condition, clearly 
from definition of y\ and 7/2 in terms of W and /^, y\ =0 and y 2 ((, 0) has an initial simple 
pole at C = Co- In this case, it is appropriate to introduce the inner variable 


where 


n v^ [1 Co 1/1 

(4.17) 

C-K 

li 

ST 

«-+. 

CH- 

(4.18) 

laOM) = 

(4.19) 

ft/,- n\ 7? 2 (Co,0) 

a = 2/2 ( 1 / Co > 0) = 

(4.20) 

fi 1 

Co/c(Co, 0 ) + Co 2 /cc(Co,0) 

(4.21) 


On substituting (4.18) and (4.19) into (2.18) and (2.26), it is clear that for y = 0(1) and 
t << 1, with 0{t ) error, we can write 


-yY{ + y\Y{^-Y 2 


(4.22) 
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Yi[r t -Y 1 ] + Y 2 [l+Y{} = 0 (4.23) 

where the prime denotes derivative with respect to 77 . Now to match to the initial condition, 
it is clear that as 77 — ► 00 , 

F 2 (r 7 ,t) ~ - (4.24) 


Y\ 


1 

V 

1 

V 


(4.25) 


Since the boundary conditions and the leading order equations do not involve t explicitly, 
to the leading order for t << 1 , each of Y\ and Y 2 will be purely a function of 7 ]. We can 
write the solution down implicitly as 

•V/2 - 1] 




, 3/2 


Y 2 = x 1 / 2 eV 2 


, , 1 f x [e*'/ 2 

f 1 = X “ 1/2 (i ~ e )- 2 i dx '~ 


1 ] 


,3/2 


(4.26) 

(4.27) 

(4.28) 


Since 7](x) i s an analytic function and 7] x ^ 0 in the finite x plane except at the origin, 
it follows from formulae (4.27) and (4.28) that the only singularities of Y 2 and Y\ are at 
X = 0, which corresponds to rj = 00 , and at x — 00 y where as Re x —cxd, 77 takes a 
finite value 

[e^ 2 - 1] 


rjo = 


-U 


d X - 


X' 


3/2 


(4.29) 


There are two possible values for i]o , depending on the choice of the branch of 1 2 . From 
numerical integration, we find 

770 = ±1.2533?'. (4.30) 

These two points are the only singularities of F(0 and F(^) in the finite 77 plane. Note 
that if Re X — * 00 , then 77 — > 00 and so at 77 = 00 , there are other possible behaviors, 
besides (4.24) and (4.25). However, this behavior occurs in the other branch sheet of the 
Riemann surface generated by the branch points 770 (given by (4.29)). To deduce the nature 
of the singularity of Y\ and Y 2 near 77 = 770 , we notice that 


X 


2 In (77 — 770 ) + 3ln ln(rj — 770 ) + 


so that for 77 


Vo, 


r, ~ 

Yj ~ 


Vo 

1 


(>) - ’fo) In (n - ijo ) + 


-(i) - Vo ) [2 In <’l - i/o)l •• 


(4.31) 

(4.32) 

(4.33) 
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These singularities at specific r/o correspond to moving singularities in the ( plane at 
Co — \Zab(orjot, as can deduced from (4.17). These singularity locations can be expected to 
be accurate only for t << 1, because of the limitation of our analysis. 

4c. One-half singularity birth at infinity 

We now assume condition (c) holds. We will show that this results in the birth of a 
one-half singularity at infinity. 

From given conditions for case (c), we get 



yi(C,0) ~ 6, C” 

(4.34) 


y 2 (C,0) ~ b 2 C m 

(4.35) 

for some constants b\ and b 2 

with n > 1 and m > 0. In this case 

, we find that the regular 

perturbation expansion in powers of t break down for ( so large 
larger. Hence, it is appropriate to introduce an inner variable 

that ( t 1 / n = 0(1) or 


v = C* 1/n • 

(4.30) 

Then upon substituting 

yi(C(*7i0»0 = \ r i(y,f) 

(4,37) 


y2(C(»M)>0 = tm/n 

(4.38) 

into equations (2.18) and (2.26) with rj = 0(1) and t « 1, 

we obtain the following 

relations (with 0(t ) error): 

-Y l + lY{-r ) Y l Y{ = 0 
n 

(4.39) 



(4.40) 

where the prime denotes derivative with respect to rj keeping t fixed. The initial conditions 

(4.34) and (4.35), translate to the following equation for j/ — > 0 



Y,( v ) - M n 

(4.41) 


Y 2 (rj) ~ b 2 rj- m . 

(4.42) 

The solution to (4.39) and (4.40) that matches with the asymptotic conditions (4.41), 

(4.42) is given implicitly as 

/y \ 

/ r l \ -Y, 

y = y- e 

\o l / 

(4.43) 
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(4.44) 


Y 2 = b? /n b 2 Y- m/n [l-nY 1 }. 

From (4.43), it is clear that r) is a regular function of Y\ everywhere in the finite complex 
Y\ plane except at the origin. Further, it is clear that Jyy has a simple zero at 1 1 = Y 
Thus, it is clear that the only singularity(ies) of Y i in the finite jj domain (other than 
tj — 0) is located where 

1 jn 

e -1 / n . (4.45) 

Note that this corresponds to n distinct locations corresponding to n distinct branches in 
(4.45). The behavior of Y\ and Y 2 near such an i]o is clearly 

Y\ ~ - - \f — 2n , / n-1 6] /n e 1 / n (r/ - t? 0 ) 1/2 + - (4-46) 

n v 

From (4.44), it follows that 

Y 2 ~ 6” ,/ ”6 2 n m / n+, y-2n 1 /— + - . (4.47) 



Thus, in the case when m > 0 and n > 1, there is instantaneous generation of n 
1/2 singularity(ies) at oo that move into the finite ( plane. The behavior of y\ and y 2 at 
each singularity is clearly seen to be consistent with (3.13) and (3.14). We find that for 
t « 1, 

C s (t)=f- 1 /"r 7o (4-48) 

Ai(t) = - - (4.49) 

Tit 

A 2 (t) = (4.50) 


B 2 (t) = 


& 6/ l/ "6 2 n m/n+1 \/-2n 1 /"- 1 6l /n e 1 / n . 


(4.51) 


5. NUMERICAL CALCULATIONS AND RESULTS 

So far, our numerical calculations liave been limited to initial conditions of the special 

form 

/(C,0) =0 (5-1) 

W{ C,0) = -eC- ( 5 - 2 ) 

This corresponds to a sinusoidal perturbation in the vortex sheet strength at an initially 
flat interface. It is likely that the numerical procedure described below can be generalized 
to other cases; however, our primary motivation in this calculation was to see if indeed the 
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approach of a one-half singularity as (3.13) and (3.14) towards |(| = 1 corresponds to a 
spike at the physical interface and illustrate through a simple example how the unphysical 
domain calculation allows us to extract all the relevant information about such a singularity. 
We chose this special initial condition for simplicity. Further, we only studied the details 
for e = 0.1, although we checked to see that similar qualitative features appeared for other 
e. From the analytical evidence presented in Section 4c, one can expect only one singularity 
forming at oo, which according to (4.37) and (4.46) will start moving down the negative 
real £ axis towards — — 1. This is confirmed by the numerical results presented here. 

We use the unphysical equations in the form (2.29), (2.32) and (2.33). The variables 
which are advanced in time are the the 6 N complex point values of ((£,t), ^ and ^ at 
N uniformly spaced out points on each of two circles in the £ -plane centered at the origin 
of fixed radii p 0 and p u where 1 < p\ < po- These set of 2 N points in the £ plane 

will henceforth be referred to as collocation points. N, chosen between 64 and 256 in all 
the calculations presented, is taken to be a power of 2 to allow convenience of fast Fourier 
Transform. 

We now describe the step by step procedure to determine the right hand side of (2.29), 
(2.32) and (2.33) at the collocation points at any time t once ((C,0> ^ ancl ^ are knnwn 
at these points. This lets us advance these variables in time, using a standard ordinary 

differential equations solver. 

Step 1: 

From the known point values of £ we use a fast Fourier Transform in Arg f to compute 
the derivative Q at the collocation points. We then evaluate and y2(C(C>*)>0> 

using the relation (2.27) and (2.28) at each collocation point. 

Step 2: 

Assuming that the image £(£,t) of |£| — p\ is outside the unit £ circle and that each 
of t/i and j/2 are analytic functions of C on and inside this curve (assumptions checked 
aposteriori), we obtain 


yi«o,o = 2 hi 

r y.(C«,0. <)<««.*) 

ki=p, C(0 - Co 

(5.3) 

92(c "’ <) = 5 hj 

f w««, <).*)<«({, 0 

I«i=p, C(0 - Co 

(5.4) 


( yi(C«,<).‘)C«K,0 

i«i=pi [C(0 - Co] 2 

(5.5) 


We compute y\(e iv ,t) and y 2 (e il ',i) for N uniformly spaced out u in the real interval 
(0,2tt). The integrals in (5.3)-(5.5) were implemented by using a real angular variable v in 
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the representation £ = p\ e lu and then using a trapezoidal rule that uses the N uniformly 
spaced points in u over the interval (0, 2 tt). 

Step 3: 

Using the computation in Step 2, Re{y\ y |] and ^|j/i| 2 are computed at N uniformly 
spaced out points on |£| = 1 circle which is then used to compute cj and Cj (for 0 < 

j < N j 2) in the Fourier representations 

OO 

Re[y\{e"' ,t)y2( etl/ >*)] = X! el}U ( 5 - 6 ) 

J — — QO 


1 OO 

-jl!/.(e , ',<)| 2 = £ (5.7) 

j~ — OO 

This is then used to calculate f) anc ^ lu (C(£i0>0 ^ ie collocation points by 

using the representations 

OO 

*3«, 0 = -«oC- 2 (5.8) 

J= 1 

OO 

n c CC,<) = 2 YtihC*-' (5.9) 

J= 1 


that follow from (2.11), (2.15) and (2.22). 

Step 4: 

Using (5.3)-(5.5) and taking £o = ^ ^ for each of the collocation points, we compute 
R 2 and i?i c by using (2.20) and (2.21). 

Once steps 1 through 4 are implemented, w’e are in a position to evaluate the right 
hand side of (2.29), (2.32) and (2.33) at the collocation points. This allows us to advance 
each of and in time by using a standard ordinary differential equations solver. 

We monitored the computed values of |C(£,t)| th e collocation points to make sure that 
they were outside the unit circle, or otherwise (5.3)-(5.5) will not be valid. Further, when 
\(\ < 1.04, we discontinued calculations since accurate evaluation of the integrals in (5.3)- 
(5.5) required a value of N larger than 256, the largest value we allowed in our calculation. 

As a biproduct of the above calculation, we get the physical interface location since 
the knowledge of j/2(C?0 on the unit circle gained in Step 2 allows us to reconstruct (up 
to a time dependent constant) the conformal mapping function on \(\ = 1 since 


,(e^=*(l,0 + /V^. (5.10) 

Note that z(l,t) cannot be found from the unphysical equations since y\ and ?/2 involve 
the derivative Zf. However, for any symmetric initial condition, z(\,t) is purely imaginary. 
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It is well known that the average height (in the physical 2 domain) f 0 ydx is conserved 
in time and so 

dv Im z(e ,v ,<) Re z v {t iv ,t) = dv Im z(t iv , 0) Re z„(e iv , 0) . (5.11) 

Since 2(£,0) is known initially, (5.11) can be used to compute 2 ( 1 , t) in (o.lO), which 
allows us to compute z(e iu ,t) completely. For purposes of computing the integrals in 
(5.10), (5.11), we used a trapezoidal rule using the N uniformly spaced out points on the 
unit ( circle, where y? is known. 

We also computed the following integrals using a trapezoidal rule and the values of 
the functions at the collocation points 



1 

27TZ 



(5.12) 



Here is obtained at the collocation points through fast Fourier transforming ( as a 
function of Arg £ on |£| = po and j^| = p\. 

It was observed that for t small enough, all the integrals (5.12)-(5.18) were zero up 
to numerical accuracy (five digits in the worst case presented here), as must be the case if 
there is no singularity of the integrands within the annular region in the £ plane between 
|£| = po and |£| = p,. Then, after some time, which was found to depend on the choice 
of po, the value of N 0 changed dramatically over a short period of time before settling 
down to a value of 1 (up to five digits). The transient time when No was significantly 
different from 0 or 1 became progressively shorter as N was increased. We interpreted this 
result to mean that around this time, a zero £ 0 (t) of crossed from |£| > po into the 
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region |£j < po . The dependence of the results on N is not expected as the numerical 

quadrature is not adequate during the short period of time when £ 0 (/) is very close to 
|£| = />o- Note however, this inaccuracy has no impact on the actual solutions to (2.29), 
(2.32) and (2.33) since it was found that the right-hand sides of (2.29), (2.32) and (2.33) 
were well behaved when this happened. This is not unexpected since there is no division by 
in implementing steps 1 through 4, as necessary to calculate the integrals (5.12)-(5.18). 
Figures 3 and 4 show the images of 256 uniformly spaced out points on |£| = po =5, at 
t = 1.4 and t = 1.6, just before and after No changed its value from 0 to 1. At / = 1.4, 
the zero of Q is at £o = —5.295 and the corresponding ( 3 (singularity location in the ( 
plane) is at -2.061. The apparent cusp in Fig. 3 is because is close to 0 at £ = —5 due 
to the proximity of ^ 0 - At a slightly later time when £o is exactly -5, we get a true cusp. 
The image of |£| = 5 in the ( plane at t = 1.6 clearly shows that the curve is not simple 
and that it intersects itself. We clearly see that the failure of one to one correspondence 
of £ and ( on this curve. Relating the geometric nature of the observed image in Fig. 4 
to general mapping properties of an analytic function, it is clear that there is a zero of 
in \(\ < po? as is consistent with the computed value of Nq = 1 at this time. We also 

monitored the image of [£| = p\ and it was found that throughout our calculations, its 
image consisted of a simple curve completely outside \(\ = 1. Figure 5 shows the image of 
256 points on j£| = p\ — 1.5 under £(£,<) at t = 2.4, Another property that was apparent 
in our calculations is that the image under £(£,/) had the effect of moving points both 
radially outwards and tangentially away from ( = 1, which corresponds to the bubble part 
of the interface, while points tended to move both radially inwards and slide tangentially 
towards £ = — 1, which corresponds to the upward moving spike (recall gravity is upwards) 
in the physical domain. This tendency is clear in Figures 3-5. There are a lot more points 
near the negative real ( axis than the positive £ axis, even when all the points were initially 
uniformly spaced out on a circle centered about the origin. Furthermore, the points near 
£ = — 1 have moved closer while image points near £ = 1 are further away than they were 
initially (Recalling £(£, 0) = £). Since the image £(£, t) defines the motion of characteristic, 
it is clear that information from the unphysical plane tends to flow towards a spike that 
will make the spike features sensitive to small changes in specified initial conditions in the 
physical domain |£| < 1. Near the bubble however, the opposite will be true. Further 
implications of this property will be discussed in Section 7. 

From the calculus of residues, it follows that for £(£,<), analytic in £, when a zero 
of Q is in the annular region the computed £ 0 in (5.13) is indeed the location £o(t) of 
the zero of Q. Further, the value in (5.14) will be the corresponding image ( 9 (t) in the 
( plane. The value zq in (5.15) is clearly C^(^o)- Further, the computed A\ y A 2 and B 2 
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in (5.16)-(5.18) has to correspond to to j/i(£o>0> J/i*(£o ,<) and 2 / 2 ( 6 ), 0- Whenever the 
calculated No = 1, the relations (3.15)-(3.17) are used to compute Aj, A 2 and B 2 in the 
expansion (3.13)-(3.14). The calculation had to be discontinued when any of the following 
conditions occurred 

(i) The image of any of the collocation points on for |£| = p\ was in |(| < 1.04 since 

this made it difficult to get a sufficiently accurate result in the quadratures (5.3)-(5.5) for 
N not exceeding 256. 

(ii) The image of any collocation points on |£| = po under the computed £(£ ,t ) came inside 
the unit ( circle as this made the evaluation in (5.8) and (5.9) sensitive to the small errors 
in Cj and cj. 

(iii) Whenever [£o| computed in (5.12) was smaller that p \ , as otherwise we could not 
ensure that the singularity £ a (t) was outside the image of |£| = p\ in the ( plane, as 
necessary for (5.3)-(5.5) to be valid. 

Since we were interested in tracking the singularity, the choice of po and pi in our 
calculations was dictated by the necessity that the calculations could be carried out over 
a significant interval of time for which No was nonzero and none of the conditions (i)-(iii) 
resulted. For all the calculations reported in the table below, we used p\ = 1.5, though we 
changed p\ between 1.25 and 2 to ascertain that the calculations of the physical interface 
did not depend on the choice of p\ as would be the case if the image of |£| = p\ in the 
( plane did not contain any singularity of yi((,t) and y 2 (C> 0 - The interfacial location 
also agreed with a standard physical domain code based on satisfying ( 2 . 8 ) and ( 2 . 12 ) 
at uniformly spaced out points on the unit £ circle, with / and W having a truncated 
representation of (2.2) and (2.3). This helped us further confirm that there was indeed no 
singularity of y\ and y 2 for |£| < |£ a |. We also ensured that the values of the integrals in 
(5.12)-(5.18) at any time t were independent of p\ when No had settled to the value 1. The 
value of po was varied from a maximum of 50 to a minimum of 3.0 in order to track the 
singularity at various stages of its motion. For a large value of po , the value of No settled 
to a value of 1 at a relatively small t. However, in such cases the computation could not 
be carried out for a long time because of the limitation of condition (ii) cited above. In 
that case, one or more points on the secondary lobe (as in Fig. 4) of the image of |£| = po 
tended to move in into |(| < 1 before long. For smaller values of p 0 in the range we 

tried, computation could be carried out for a significantly longer time; however for such 
cases, singularity tracking was not possible for earlier times since £o{t) was then outside 
|£| = p 0 . However, smaller p 0 allowed us to track ( s in the later stages when ( s was getting 
fairly close to -1 when the corresponding interface shape near £ = — 1 showed a developing 
spike. We also noted that for two different po, in the overlapping time interval when the 
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computed No had settled to a value of 1 , there was agreement in the computations (5.12)- 
(5.18). This provided an additional check on the code and on our assumption that there 
was indeed only one singularity in this case. 

Table 1 lists the various quantities of interest as a function of time. The results have 
been checked by appropriately doubling N until there was no variation in the results. 
With the limited number of N < 256, we had difficulty ensuring desirable accuracy for 
t greater than 2.4. The physical interface at t = 2.4 is shown in Fig. 6 . The image 
of ( = —1 corresponds to the spike observed at x — n. The curve marked 1 in figures 
7, 8 , 9, 10 show i4i, A 2 and Bz as a function of t, where we used (3.15)-(3.17) and 
the computed quantities in (5.12)-(5.18). The curve marked 2 is the analytic prediction 
(4.48)-(4.51), where b\ = — e = — 0 . 1 , 62 = 1, ti — 1 and m = 0. Since the theory in 
Section 4 requires t << 1 , we get surprisingly good agreement even for t not all that 
small. Figure 11 shows Im z(l,t) and 7m z(— l,t), the vertical location of the bubble tip 
and spike respectively. Though, the spike does not appear very well developed in Fig. 6 , 
Fig. 11 shows that that the long time asymptotic range has been reached where the spike 
accelerates upwards (the direction of gravity) with the acceleration approaching free fall. 
Fig. 12 shows the product Re[A 2 B%] as a function of time. In this case, since each of Az 
and i ?2 are real, Re [A 2 B 2 ] = AzBz- This quantity appears to reach a minimum and then 
increase. Thus M\{t ), as defined in (6.32), goes through a maximum and then decreases. 
We take this as an indication that M\(t) does not blow up in finite time. We cannot rule 
out the possibility that the observed trend in M\{t) reverses at at an even longer time; but 
given that the spike has reached its asymptotic acceleration towards the the end of our 
calculations, we do not expect this to happen. In the following section, analytical evidence 
will be presented to show that under the assumption that f* M\(t)dt does not blow up in 
finite time, an isolated one-half singularity of the type (3.13) and (3.14) cannot reach the 
physical domain |£| = 1 in finite time. 

6. EVENTUAL FATE OF A ONE-HALF SINGULARITY 


We now want to understand the fate of an isolated singularity of the type given by 
(3.13) and (3.14), once it is already close to the unit circle in the ( plane. 

Note that (2.21), (2.22) and (3.18) imply that 


^ = -h - y 2 (^l^,0 yi(C*(0>0- 


Now, from (2.25), (3.13) and (3.14), it follows that 


h(Cs,t) 


jl r d „ 

2 ^ Jo 


e' v + C» 

t iv - (s 


Q(v,t) 


( 6 . 1 ) 


( 6 . 2 ) 
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where 

Q(M) = He [»,*(€“, t),t) ,t) ] . (6.3) 

It is convenient to define 

Ql(M) = Re [x;(t)B 2(0 (=■'■' -C.(‘))' /2 ] (6.4) 

<? 2 (M) = He [Af(t) B 3 (0 («“ - 0,(0) ] (6.5) 

Qs(»,0 = Re [AJ(t) B 2 (()] |e ,v - C.WI- ( 6 - 6 ) 

If we denote 

0 = R <?“' , (6.7) 

then for R — > 1 " , the behavior of Q(v,t) in the vicinity of i / 0 is given by 

Q(v,t) = Qi(M) + Q 2 (M) + £h{v,t) + Odi' - ^o| 3/2 ) • ( 6 . 8 ) 

Now let’s define 

1 f 2n S e iv 4- C 

/ 21 ((„0 = ^ J dv e , v _ c ” Q .(»,0 (6-9) 

1 Z’ 2 * re" 7 4- f 

/««.,<) = 5 - / 77 _ , « 2 (".') ( 6 - 10 ) 

■^ 7r Jo L e ss . 

1 /“ 27r re’ 17 4- f 

MC .,0 = ~ dv S I T *(> 60 . ( 6 . 11 ) 

- 7r Jo L e s* . 

By changing the integration variable from i/ to — i/ in each of (6.9), (6.10) and (6.11), 
it is not difficult to establish the property that for 1 < j < 3 , 

1 e iv -r- 

&*«.,<) = -=" / >!» fc 0 2 (-> 60 . ( 6 . 12 ) 

j o e — <7 

Now note that from previous expressions (6.4), (6.5) and ( 6 . 6 ) 

<?,(-*, 0 = He [.4, (t) B-(t) (r'* -C(t))' 11 ] (6.13) 

<M-M) = He [^(0H,*(0(e'"-C(0)] (6-14) 

<?.,(-», 0 = He [A 2 (0 B 2 ’(01 k'" - C(0I- (6.15) 

Each of (6.13) and (6.14) are clearly seen to be the real part of simple analytic functions 
for £ = e'" . Therefore, it immediately follows that 

£. «.(<). o = -x,(o b 2 *(o - cm) (6.i6) 
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fo(C,(t).i) = -Mt) b M(£m - CM)- 


Now consider I 2 3(( a (t),t) ■ From (6.6) and (6.11), it follows that 

h j(C.(0,«) = Re [Al(t) B 2 (t)] j dv \t* - c;(0l 


e + C.(0 

^ iv \ r 

e crm 


(6.18) 


Using (6.7) and periodicity and symmetry of the integrand, it follows that 


I 23 — 


(■ R 2 - i) 


Re [. A*(t ) fl 2 (t)] 


(i? 2 4- 1 — 2 i? co-si/) 1 


(6.19) 


Now we want to consider the asymptotics of (6.19) in the limit R — ► 1 . To do this it 

is convenient to break the integral range in (6.19) 


l = l + l 


( 6 . 20 ) 


where 1 >> e >> (R — 1) . For the first integral in (6.20), the leading order asymptotics 
can be found by replacing R by 1. For the second integral in (6.20), the leading order 
asymptotics can be found by first simplifying the cos v by 1 2 and then replacing 

the integration variable v by u , where v = (R — 1) v . One then finds the leading 
order asymptotics of (6.19) as 




-2 Re ( A*B 2 ) 

7T 


* dv 
2 sin | v 


dv 

) (1 + V 2 ) 


1/2 ' 


( 6 . 21 ) 

Each of these integrals can be computed exactly with e dropping out completely on 
addition. We find 


l23((s(t), t) 


= -2 (JLJ1 Re (AS(<) B 2 (t)) [3 In 2 - In (R - 1) + o(l) ] . (6.22) 


Now let compute 


i R (Ui),t) = - /21 (C.(0i0 - Jm(C*(0>0 - ^23 (C.(0.0- ( 6 - 23 ) 


It is clear that 


where 


««*)■*) = h jf 


, \e iv + C. 1 ^ ^ 

dv — — — Vk(M) 

e . 


Qfl(v,t) = Q(l/,t) - Q\(v,t) - Q 2 (v,t) - Qz{v,t). 


(6.24) 
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This can also be written as 


/*«•(*)>*) 


1 f 2ir R 2 + 1 — 2iRsin{v — v 0 ) 

2tt J 0 1 + R 2 - 2Rcos(v - vo) 


Qk(M)- 


(6.26) 


As -> 1 , the above expression reduces to 




1 

7T 



dv 


Qro (M) 

2 — 2cos(v — vo) 


[1 + 0 ( 1 )] 


- f dv cat \ {v - Vo) Qr 0 {v,t) [1 + o(l) ] (6.27) 

ZTT J o Z 

where the subscript 0 in Qr refers to its evaluation with () s = <’ !</ ° . Notice that from 
the behavior (6.8), each of the integrals in (6.27) exists. 

Now note the fact that if for any complex £ s singularity, near which the asymptotic 
behavior of j/i and y 2 is given by (3.13) and (3.14), there is a singularity at in the 
neighborhood of which 


»,k,o = a\(o + AW)«-c(t ))' /2 + Aim -a + »(c - o 3/2 - 


(6.28) 


and 

tte(C .<) = B’ 2 (t) (( - C(<))' /2 + BSm - C) + <*« - C) 3/2 )-. ■ (6-29) 

Now as R — * 1~ j the point ^ comes close to £ = Cs s-nd so the last term in (6.1) 
yields 

1 . 1 . 1/2 
-y 2 (7777^) yi(C*( 0 » 0 ~ a O) b 2 (7-777 - CO) 

CO CO . (6.30) 

+ .Ai(t)i?3 (— y — co) + o(i 7 — 1) 

Taking the real part of (6.1), using (6.7), (6.16), (6.17), (6.22), (6.27) and (6.29) in the 
asymptotic limit R —> 1~ , we get 

R = M\(t) (R — 1) ln(R — 1) + M 2 (t)(R- 1) + o(i?-l) (6.31) 


where 


M 2 (t) 


M\(t) = --Re [A*(t) B 2 (t) } 

7T 


— Re [A3(0 B 2 (()] 


I f’ dv . 

7T J 0 2-2 cos (v- Vo) 


(6.32) 

(6.33) 
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The solution of (6.31) to the leading order is 

ln(R— 1) = In (i?o _ 1) ^ \ (6.33) 

If we assume that j* M\ ( t)dt does not blow up for any finite time, then it follows from 
(6.34), that there will be absence of finite time singularity. This assumption is consistent 
with the results from numerical calculations as reported in the previous section, though 
for a special initial condition. In that case, recall that M\ (i) actually appeared to decrease 
after sufficiently long time. 

7. DISCUSSIONS AND CONCLUSION 

We have analyzed some properties of the unphysical equation derived earlier by Tan- 
veer (1991a). For early times, for a specific class of analytic W((,0) and z((,0) in |C| > 1 
(except possibly at infinity), we noted how one or more singularities can form at a point in 
1(1 > 1, where 2f((,0) is zero. One-half singularities are shown to be generated, though 
we show that other singularities of a much more complicated form involving logarithms 
are possible. Further, analytical evidence suggests that one or more one-half singularity 
may be born at ( = oo, which moves to the finite ( plane instantaneously. The numerical 
computation, although for a special case, clearly shows that one-half singularities on their 
approach towards the physical domain corresponds to a continually developing spike. The 
connection of more complicated singularities with interfacial features is yet to be made. 
The numerical computation also suggests that the characteristic field in the unphysical 
plane is pointed away from the bubble and towards the spike. With an assumption that 
appears to be consistent with the numerical calculation for the special case, our analytical 
evidence suggests that an isolated one-half singularity of the type given in (3.13) and (3.14) 
cannot reach the physical domain in finite time. We have not investigated the question of 
finite time singularity in the physical domain when possibly multiple one-half singularities 
merge or other forms of singularities are present, and this remains an open question. 

It is clear that the location of formation of initial singularities depend on the specifics 
of the initial conditions ^((,0) and W((,0) in |(| > 1. Clearly, it is possible to make 

arbitrary small perturbation in the physical domain |(| < 1 that significantly alters the 
location and number of singularities that are formed in |(| > 1. Further, our numerical 
computation, though for a special case, appears to indicate that the characteristic field in 
the unphysical plane is directed away from the bubble and towards the spike. This feature 
accounts for the sensitivity of spike to initial conditions in the physical plane. Further, if 
one were to perturb the initial condition (5.1) and (5.2) by placing an additional singularity 
in |( | > 1 that is weak in the sense that its contribution to the right hand side of (2.29) 
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is small, then the characteristic field generated in this problem will be close to what has 
been computed in Section 5. Since a singularity moves along with the characteristic field, 
one can expect that in due time, this additional singularity will also approach ( = -1, 
which corresponds to the spike. 

This suggests that in a general way, it might be true that weaker singularities have 
the tendency to merge with stronger ones. If this is generally true, we speculate this as a 
possible reason why bubble competition results in a dominant bubble. The image under 
of an arc on |(| = 1 between two approaching singularities must contain the bubble 
region, since each spike possibly corresponds to a singularity. If the singularities merge 
pairwise in the ( plane, the corresponding bubble region between these two singularities 
will disappear resulting in one larger bubble. However, if multiple singularities of nearly 
equal strength results from an initial condition, this merger can be expected to t<.ixe a 
while as the characteristic field will be almost equally affected by all the singularities. The 
transient dynamics observed in experiment can be expected to depend on the location of 
initial singularities and their motion, the randomness resulting from the ill-posedness in 
determining z((,0) and TT((,0) in |C| > 1, when they are only given in the physical 

domain |(| < 1. Confirmation of this scenario must await further investigations in the 

future. 
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t 


C. 

A\ 

A-2 

B 2 

Spike 

Bubble 

0.24 

-41.27 

-15.20 

4.249 

— 1.541 

0.3626 

0.024527 

-0.023940 

0.40 

-24.38 

-8.994 

2.639 

— 1 .244 

0.4721 

0.041930 

-0.040243 

0.60 

-15.72 

-5.820 

1.868 

-1.095 

0.5843 

0.065733 

-0.061680 

0.80 

-11.28 

-4.201 

1.516 

-1.046 

0.6857 

0.092863 

-0.084978 

1.00 

-8.533 

-3.209 

1.327 

-1.048 

0.7808 

0.12467 

-0.11087 

1.20 

-6.657 

-2.539 

1.221 

-1.083 

0.8712 

0.16283 

-0.14008 

1.40 

-5.295 

-2.061 

1.161 

-1.142 

0.9555 

0.20949 

-0.17333 

1.60 

-4.268 

-1.711 

1.132 

-1.220 

1.029 

0.26745 

-0.21131 

1.80 

-3.477 

-1.454 

1.126 

-1.315 

1.085 

0.34039 

-0.25464 

2.00 

-2.860 

-1.269 

1.141 

-1.426 

1.108 

0.43316 

-0.30381 

2.20 

-2.381 

-1.141 

1.181 

-1.559 

1.080 

0.55200 

-0.35913 

2.40 

-2.013 

-1.062 

1.248 

-1.729 

0.9805 

0.70441 

-0.42066 


Table 1: Various parameters characterizing the one-half singularity. The last two 
columns contain Im and Im z(l,t), the vertical locations of the spike and 

the bubble. 
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Fig. 1: The complex ( cut unit circle. 
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Fig. 2: The physical z = x + iy plane. 


31 



Fig. 3: The image of 256 uniformly spaced out points on |£| = p 0 — 5 at t = 1.4, just 
before the crossing of £ 0 into |£| < 5. 
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Fig. 4: The image of 256 uniformly spaced out points on |£| = po — 5 at t — 1.6, a 
little after £o has crossed into |£| < 5. 
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Fig. 5: The image of 256 uniformly spaced out points on |£| = p\ 
at t = 2.4, which remains a simple curve. 


3.0 

1.5 in the ( plane 
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x/( 2 pi) 


Fig. 6: The shape of the interface at t — 2.4 when normalized by the wavelength 2n. 
A spike can be noted at = 0.5 corresponding to ( = —1. 


zeta 
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I 



'O 


0.0 





t 


Fig. 7: The trajectory of the singularity ( a (t) as a function of t is shown by the 
curve marked 1. Curve 2 correspond to the analytical prediction (4.45) using a small 
t analysis. 
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.5 1.0 1.5 2.0 2.5 

t 


s a function of t in the curve 1. Curve 2 corresponds to the 
; small t analysis. 
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c4 

0.0 0.5 1.0 1.5 2.0 2.5 

t 

Fig. 9: A 2 is shown here as a function of t in the curve 1. Curve 2 corresponds to the 
analytical prediction using small t analysis. 
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Bubble and Spike location 
- 1.0 - 0.5 0.0 0.5 



0.0 0.5 1.0 1.5 2.0 2.5 

t 


Fig. 11: Plot of Im z(\,t) and Ira z( — l,t), the vertical location of the bubble-tip 
and the spike as a function of t. 
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